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ABSTRACT 

[Secret] 

Several  statistical  topics  were  studied  concerning  detection  of  tar¬ 
gets  on  the  ocean's  surface  from  a  sate  Hite -based  radar.  First,  a  de¬ 
sired  signal-to-noise  ratio  was  obtained  by  considering  the  tradeoff 
between  false  alarms  and  detected  targets.  For  a  0.99  detection  prob¬ 
ability  and  a  10'  false  alarm  probability  the  S  N  ratio  would  be  16  db 
for  a  nonfluctuating  target.  Second,  the  decorrelation  times  of  sea 
clutter  were  calculated  and  found  to  be  so  small  that  for  side-looking 
radars  sea  clutter  samples  are  independent  from  pulse  to  pulse  but 
found  to  be  so  large  for  forward -scanning  radars  that  one  finally  ob¬ 
tains  decorrelation  from  sea  motion  instead  of  from  platform  motion. 
Third,  a  fluctuating  target  was  considered,  and  the  losses  were  calcu¬ 
lated  and  compared  with  those  of  a  nonfluctuating  target.  Roughly,  one 
needs  1  to  4  db  more  power  to  detect  a  fluctuating  target  as  opposed  to 
a  nonfluctuating  target,  Th^n,  the  ionosphere  was  considered,  and  the 
effects  of  Faraday  rotation  a.  d  the  random  phase  shift  were  calculated. 
Degradation  varies  from  Intolerable  below  900  MHz  to  negligible  at 
3  GHz.  Next,  some  of  the  aspects  of  data  processing  were  studied;  the 
optimal  integration  angle  was  found,  the  optimal  weighting  was  calcu¬ 
lated,  and  azimuthal  position  estimators  were  considered.  Surprisingly, 
optimal  weighting  gives  a  pattern  only  0.3  db  better  than  uniform 
w’eighting.  Finally,  a  wav  was  determined  of  using  an  adaptive  threshold 
and  soft  limiting  to  maintain  a  constant  false  alarm  rate. 


PROBLEM  STATUS 

fills  is  a  final  report  on  one  phase  of  the  problem;  work  on  other 
phases  continues. 
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OCEAN  SURVEILLANCE  STATISTICAL  CONSIDERATIONS 
(Unclassified  Title) 


SUMMARY  (S) 

After  the  development  of  ground-based  operational  naval  radars  in  the  1930's,  they 
were  placed  on  ships  and  later  in  aircraft.  Now  that  the  space  age  is  here,  a  logical 
evolution  should  include  placing  a  naval  radar  in  a  satellite.  This  report  discusses  some 
of  the  general  statistical  considerations  associated  with  the  detection  of  targets  on  the 
ocean's  surface  from  a  satellite-based  radar. 

In  the  first  section  the  probability  of  false  alarm,  the  probability  of  detection,  and 
the  required  signal-to-noise  ratio  for  ocean  surveillance  are  discussed.  From  calcula¬ 
tions  of  the  number  of  false  alarms  one  would  obtain  iier  orbit  it  appears  desirable  to 
have  a  very  low  false  alarm  rate,  i.e.,  10  -’  to  lO"'".  Thus,  if  a  probability  of  detection 
of  C.99  is  desired,  an  integrated  signal-to-noise  ratio  of  about  16  db  is  necessary  for 
nonfluctuating  targets. 

In  the  second  section  the  decorrelation  times  of  the  sea  clutter  due  to  the  radar’s 
platform  motion  are  calculated.  The  decorrelation  times  for  a  side-looking  radar  are  so 
small  that  one  has  independent  sea  clutter  samples  from  pulse  to  pulse.  However,  for  a 
forward-scanning  system  the  decorrelation  times  are  so  large  that  one  finally  obtains 
decorrelation  from  the  sea's  motion  as  opposed  to  the  platform’s  motion. 

In  the  third  section  a  simplified  target  model  consisting  of  a  collection  of  point  scat- 
terers  is  proposed.  Employing  this  model,  a  Monte  Carlo  method  is  used  to  calculate 
the  number  of  independent  target  returns  one  receives  from  this  target  model  as  a  func¬ 
tion  of  integration  time,  frequency,  and  aspect  angle.  Then,  by  using  Swerling’s  curves, 
the  required  signal-to-noise  ratios  to  olitain  a  probability  of  false  alarm  of  10* and  a 
probability  of  detection  of  0.9  for  fluctuating  targets  is  found  as  a  function  of  frequency 
and  integration  time. 

In  the  fourth  section  the  average  daily  worst  case  ionosphere  is  assumed.  From 
this  model,  the  degradations  in  the  signal-to-noise  ratios  and  degradations  in  the  signal- 
to-clutter  ratios  are  calculated.  .All  frequencies  below  900  MHi  have  intolerable  degra¬ 
dations  and  so  does  900  Me  above  200  naut  mi  altitude.  The  degradation  disappears 
swiftly  as  the  frequency  is  increased;  so  that,  there  is  essentially  no  loss  at  3  GHz. 

Also,  the  ionosphere  imparts  a  random  phase  shift  on  the  signal  wliich  restricts  the  co¬ 
herent  integration  time  to  the  time  required  for  the  ionosphere  to  shift  the  signal  0.8 
radian. 

In  the  fifth  section  the  optimal  integration  angle  is  found  to  vary  between  0.92  and 
0.88  of  the  3-db  beamwidth,  and  the  associated  scanning  loss  varies  between  1.7  and 
1.4  db  — the  variation  being  a  function  of  the  number  of  pulses  integrated.  Then,  the  op¬ 
timal  weighting  pattern  is  calculated  and  the  surprising  result  is  found  that  the  optimal 
pattern  is  only  0.3  db  better  than  the  uniform  weighting.  Finally,  several  methods  are 
presented  for  estimating  the  azimuthal  position.  These  methods  give  results  which  are 
about  20''f  larger  than  the  smallest  theoretical  accuracies  obtainable. 

in  the  final  section  a  method  of  using  an  adaptive  thresiiold  and  soft-limiting  to  ob¬ 
tain  a  constant  false  alarm  rate  is  presented.  The  soft-limiting  loss  is  shown  to  be  very 
small  if  the  limiting  level  is  about  the  3a  poii.i  of  the  noise  distribution. 
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SIGNAL-TO-NOISE  RATIO  (S) 

Single-Look  Radar  (S) 

Before  one  can  determine  the  minimum  acceptable  sigral-to-noise  ratio,  one  must 
first  decide  where  one  should  operate  on  the  detection  curve  tor  a  given  .s  n  ratio;  that 
is,  one  must  find  the  probability  of  false  alarm  ,  and  the  probability  of  detection 
For  example,  with  a  13 -db  >  \  ratio,  one  can  achieve  a  =  10"^'  and  /',  =  0.9;  how¬ 
ever,  one  could  also  achieve  a  ,  =  lO'-  and  a  corresponding  ft.  =  0.8.  To  decide 
which  of  these  two  conditions  is  better,  a  Bayesian  approach  can  be  adopted;  that  is,  a 
cost  function  can  be  generated. 

A  particular  cost  function  can  be  based  on  the  argument  that  if  one  is  trying  to  track 
the  ships  on  the  ocean,  then  either  missing  a  target  or  detecting  a  false  one  will  cause 
equal  difficulty  with  the  tracking  procedure.  The  optimal  operating  point  can  be 

considered  to  be  the  point  that  maximizes  yi/‘,  ..  r,.  s  Vi  for  a  given  v  ratio,  where 

N  I  No  ('f  I  .1  r  LM  I  s  (1 1  t  C  I  »’  tl  —  \  II  f  r  ,1  !  ii  I  .1  r  Iti  s 

}' ,  \'  n  if  I  ;i  r  c  t  I  ^  •  —  /  *.•  .  '  ^  .  of  <1  f  I  I  s  I  I  >  n  ' 

This  cost  function  is  usually  called  the  Ideal  Observer.  Middleton  (1)  discusses  this 
function  further. 

To  calculate  the  optimal  operating  point,  particular  radar  configurations  must  be 
assumed  so  that  the  number  of  targets  and  number  of  decisions  can  be  calculated.  The 
radar  configurations  assumed  are  given  in  Table  1. 


Table  1 

Radar  Configurations  Assumed 


Configu¬ 

ration 

Swath 
Width 
(naut  mi) 

Resolution  (ft) 

Range 

Azimuth 

1 

800 

200 

200 

2 

1000 

50 

1000 

3 

1000 

50 

12,000 

For  configurations  1  and  2  the  number  of  decisions  is  1.5x10'"  and  the  number  of 
targets  is  approximately  600  in  a  90-niinute  interval  —  a  typical  orbital  time.  (The  num¬ 
ber  of  decisions  is  equal  to  the  area  pul  under  surveillance  per  orbit  divided  by  the  area 
of  tlie  resolution  cell.)  Consequently,  to  find  the  optimal  operating  point  for  each  >'  ^ 
ratio,  Table  2  was  generated.  This  table  was  generated  for  a  nonfluctuating  target  from 
curves  taken  from  Ref.  2. 

From  Table  2,  one  can  easily  extract  the  optimal  operating  point  (the  point  that 
maximizes  the  cost  function)  for  each  >'  ''  ratio.  The  operating  points  for  each  >  \  ratio 
are  presented  in  Table  3. 
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Table  2 

Evaluation  of  the  Cost  Function  (Number  of  Targets  Detected  Minus  the  Number 
o[  False  Alarms)  for  Configurations  1  and  2  of  Table  1,  for  which  the  Number 
of  Targets  is  600  and  the  Number  of  Decisions  is  1.5x  lo'” 


Signal-to- 
Noise  Ratio 
(db) 

Probability 
of  a  False 
Alarm 

No.  of 
Targets 
Detected 

No.  of 
False 
Alarms 

'1 

Cost 

Function 

. 

13 

10  ■ 

480 

1500 

-1020 

10-' 

420 

150 

270 

10-'' 

330 

15 

315 

10-'“ 

240 

2 

238 

10-" 

162 

0 

162 

14 

10-^ 

570 

1500 

-930 

10-“ 

540 

150 

390 

10-'' 

480 

15 

465 

10-'“ 

420 

2 

418 

10- '  ' 

33C 

0 

330 

15 

10-- 

594 

1500 

-906 

.0- 

588 

150 

438 

10- 

570 

15 

555 

1 

10- 

540 

2 

538 

1 

10-  " 

480 

0 

480 

16 

10-- 

!  600 

1500 

-900 

10  ■“ 

599 

150 

449 

10-“ 

597 

15 

582 

1  10''" 

1 

1  594 

2 

592 

i 

10- " 

1  582 

0 

582 

i 

10-' 

1 

1  600 

1500 

-900 

1 

i 

10-“ 

j  600 

150 

I  i 

1  450 

1 

i 

10-'-’ 

i  600 

1 

15 

585 

1 

1 

10-'“ 

1  600 

2 

598 

j 

10-  " 

j  599 

0 

599 
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Table  3 

Optimal  Operating  Points  as  Taken  from  Table  2 


Ratio 

(db) 

" 

No.  of 
False 
Alarms 

No.  of 
Targets 
Detected 

13 

lo-'* 

15 

.55 

330 

14 

10  ■' 

15 

.8 

480 

15 

10-" 

15 

.95 

570 

16 

10-'" 

2 

.99 

594 

17 

10 

0 

.999 

599 

Obviously,  5  \  ratios  of  13  and  14  db  are  not  acceptable  because  of  the  low  proba¬ 
bilities  of  detection.  On  the  other  hand,  any  >  \  ratio  above  15  db  is  acceptable.  Since 
there  is  some  gain  in  going  from  15  to  16  db  but  very  little  gain  in  going  from  16  to  17  db, 
it  appears  that  16  db  is  a  sufficiently  large  >  ^  ratio. 

To  show  that  this  result  is  insensitive  to  the  particular  system  considered,  let  us 
consider  configuration  3,  where  the  azimuthal  resolution  is  12,000  ft  and  the  number  of 
decisions  has  been  reduced  to  1.2  xlO".  The  cost  function  for  this  system  is  given  in 
Table  4. 

Again  looking  at  the  optimal  operating  points,  one  comes  to  the  previous  conclusion; 
i.e.,  any  >  \  ratio  greater  than  15  db  is  acceptable.  Thus,  us  long  as  one  has  large 
swath  widths  and  small  range  resolution  cells,  an  n  ratio  of  about  16  db  will  be 
required. 


Multilook  Radars  (S) 

In  the  previous  paragraphs,  the  assumption  of  a  single-look  radar  has  been  made; 
however,  in  this  subsection  the  consequences  of  the  radar  being  able  to  look  at  each 
resolution  cell  several  times  will  be  investigated.  For  simplicity,  we  initially  consider 
a  two-look  radar,  with  the  time  between  looks  being  90  seconds.  (For  a  target  to  be  de¬ 
tected,  it  must  be  detected  on  both  looks.)  Now,  if  targets  of  interest  liave  a  maximum 
velocity  of  30  kiiots,  they  can  travel  4560  ft  in  90  seconds.  Consequently,  on  tlic  second 
look,  the  target  can  appear  in  any  of  the  following  number  of  resolution  cells  (a  cell  is 
200  ft  by  200  ft)  which  surround  the  initial  observation  of  the  target: 


If  V  is  the  number  of  false  alarms  per  look,  the  probability  that  a  false  alarm  in  the  first 
look  will  be  matched  with  any  of  tlie  1630  cells  surrounding  every  false  alarm  for  the 
second  look  is  approximately 


I', ii)  \ 
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Table  4 

Evaluation  of  the  Cost  Function  for  Configuration  3  of  Table  1,  for  which 
the  Number  of  Targets  is  600  and  the  Number  of  Decisions  is  1.2x10'* 


14  i  10-'  I  570  120  450 


10- ”  I  540  12  I  538 

I  i  ' 

I  lO*'-"  j  480  .  1  I  479 

;  lO-'-’  ,  420  I  0  i  420 

I  10- "  '  330  I  0  i  330 

I  .  ; 

15  i  10-'  !  594  i  120  !  474 

j  10-‘*  I  588  j  12  1  576 


( 

1 

10-'* 

j  570 

1 

569 

I 

10-'*’ 

1  540 

0 

540 

! 

10- " 

480 

0 

480 

16  1 

10-' 

600 

1 

120 

480 

i 

1 

10-*' 

599 

12 

587 

i 

lO-" 

1  597 

1 

596 

! 

t 

10- IV 

!  594 

0 

594 

1 

10- ' ' 

j  582 

0 

582 

17  j 

10-' 

!  600 

120 

480 

' 

10-” 

1  600 

12 

588 

10-'* 

600 

1 

599 

1 

10-  p) 

600 

0 

600 

_ L 

10- " 

599 

J _ 

0 

599 
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Finally,  the  total  number  of  false  tracks  is  approximately 

1  f>  ID  V- 
1.3^  10'-' 

If  one  can  tolerate  ten  false  alarms  per  orbit,  the  number  of  false  alarms  i^er  look, 

V,  is  9560,  and  the  probability  of  false  alarm  per  look  is 

o'.iiO 

l\  -  -  0  .  |:i-  . 

“  I .  ■  I 

If  one  requires  au  overall  probability  of  detection  t)f  .98,  the  probability  of  detection  per 
look  is  .99.  For  /’, ,,  =  0.6 >  10"''  and  =  .99  the  \  ratio  required  is  about  14.2  db. 

One  can  perform  a  similar  analysis  allowing  one  or  100  false  alarms  instead  of  ten. 

Then,  the  corresponding  \  ratios  are  14.5  and  14.0  db  respectively.  Thus,  by  using 
two  looks  instead  of  one  look,  one  can  pick  up  only  about  1.5  db.  The  gain  is  smaller  tluii' 
one  would  expect  because 

1.  A  small  change  in  tlie  \  ratio  causes  a  large  percentage  change  in  r,  . 

2.  The  targets  are  moving;  thus,  one  has  a  multitude  of  ways  to  match  up  tlie  possi¬ 
ble  targets. 

By  a  similar  line  of  reasoning,  the  number  of  false  tracks  m  i  (for  an  I'-look  radar, 
with  each  look  separated  by  90  seconds)  is  approximately 


I  I  3  .  1..'  J 

If  N i  I  I  and  '•/  1,  \  =  43,900  false  alarms  jjer  look;  and  the  corresponding  >  3  ratio 
for  a  /■  =  0.993  is  14.3  db.  Thus,  a  three-look  radar  is  only  0.3  db  better  ti.an  a  two- 
look  radar.  In  fact,  the  performance  will  improve  very  little  as  the  number  of  looks  is 
increased;  for  instance,  a  ten-look  radar  is  only  0.3  db  better  than  a  three-look  radar. 


SEA  CLUTTER  DECORRELATION  TIME  (S) 

To  determine  the  gain  associated  with  integrating  in  tiie  presence  of  sea  clutter,  a 
correlation  criterion  is  defined  as  follows:  Since  the  urea  of  the  resolution  cell  is  very 
large,  the  returned  signal  is  the  resultant  oi  scattering  from  many  rc..tterers  within  a 
single  resolution  cell,  then,  two  samples  from  a  given  resolution  cell  will  be  uncorrelated 
if  the  maximum  change  in  the  one-way  range  separation  of  several  scatterers  exceeds 
half  a  wavelength.  In  the  following  development,  the  sea  clutter  is  assumed  to  be  uni¬ 
formly  distributed  in  a  range  cell  and  the  decorrelation  is  assumed  to  be  entirely  due  to 
the  motion  of  the  radar  platform;  that  is,  the  sea  is  assumed  to  be  motionless.  At  time 
I)  (Fig.  1)  the  ranges  to  the  two  scatterers  are  equal;  at  time  ■  /'  the  r.inges  are 

given  by 

/.■;  -  ^  J/.M  .  .  (1) 

and 
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Fig.  1  -  Squinlfd  radar  conl'igur.ilion 


where 


and 


/O  i«tH 

(3) 


I  ,  OH  D)  >  OH  (  ■  .  .  OH  D) 


in  which  ,  is  the  angle  between  the  sight  vector  and  the  velocity  vector  of  the  radar  plat 
form,  D  is  the  depression  angle,  is  the  squint  angle,  and  :  is  the  beamwidtli.  Sub¬ 
tracting  Eq.  (2)  from  Eq.  (1)  yields 


-  It  :  titii  ( .  ..H 


or 


(  A’l  <  UJ  (  /i,  A’  J 
L’/.M  I  os  ■  -  .  *.s  .  ,1 


One  can  now  make  the  substitutions 


and 


t  It  ,  nt 


yielding 


A’,  It,  \  i  . 


_ I _ 

l!  (  .  os  ‘  -  .■  OH  .  ,n 
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Finally,  letting  v  -  VT.  (where  i'  is  the  velocity  of  the  radar  platform)  and  solving  for  7.,  , 
one  obtains 


Note  that  c  and  ■  must  have  the  same  distance  dimension.  Equation  (4)  was  programed 
for  the  CDC  3800  computer  at  NRL,  and  T,  was  calculated  as  a  funclion  of  frequency, 
beamwidth,  squint  angle,  and  depression  angle.  Since  the  amount  of  output  data  is  enor¬ 
mous,  Only  two  special  cases  will  be  presented:  the  side-looking  radar  and  the  forward¬ 
scanning  radar. 


Side-Looking  Radar  (S) 

When  the  radar  is  a  side-looking  radar,  the  squint  angle  equals  90  degrees.  Then 
Eq.  (3)  implies  .  =  90  degrees.  If  one  uses  the  small  angle  approximation  for  ■  ,  Eq.  (4) 
reduces  to 


/  -  . 

Cl'. 

Using  I  =  25,000  ft  sec,  / .  was  calculated  for  various  beamwidths  and  frequencies  se¬ 
lected  from  the  available  radar  bands  (Table  5).  Reviewing  the  table,  one  notices  that, 
except  for  one  or  two  instances,  all  the  returning  pulses  can  be  considered  independent. 
Consequently,  the  integration  gain  for  the  sea  clutter  is  the  same  as  for  the  thermal 
noise. 


Table  5 

Decorrelation  Time  for  A  Side-Looking  Radar 


Frequency 

I  /.  (sec)  for  Various  Beamwidths 

(MHz) 

!  0.5 

J.  .  .. 

1.0 

2.0 

140 

0.0322 

0.0161 

r-^ToTT 

0.0081 

220 

0.0205 

0.0102 

0.0068  1 

0.0051 

440 

0.0102 

0.0051 

0.0034  j 

0.0028 

900 

0.0050 

0.0025 

0.0017  ' 

0.0013 

1300 

] 

0.0035 

0.0017 

0.0012 

0.0009 

2900 

0.0016 

0.0008 

0.0005  i 

0.0004 

5250 

0.0009 

0.0004 

0.0003  ! 

0.0002 

8500 

0.0005 

0.0003  1  0.0002  i 

1  _  ! 

0.0001 

Forward-Scanning  Radar  (S) 

When  the  radar  is  a  foi-ward-scamiing  radar,  the  squint  .ipglc  equals  0  degree. 
Then  Eq.  (3)  shows  that  .'v  If  one  uses  the  small  angle  approximation  to  ,  Eq.  (4) 
reduce.s  to 
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r 

Using  P  =  20  degrees  and  v  =  25,000  ft  sec,  7 .  was  calculated  for  the  beamwidths  and 
frequencies  of  interest  (Table  6).  The  decorrelation  times  in  Table  G  are  quite  large.  In 
fact,  the  decorrelation  due  to  the  motion  of  the  sea  is  less  than  most  of  the  times  in  Ta¬ 
ble  6.  This  decorrelation  time  is  given  by 


where  i  is  the  sea  particle  velocity.  Using  Table  7  (3),  for  sea  state  4,  Eq.  (5)  reduces 
to 


where  ■  is  the  wavelength  in  meters  (  =  30  frequency  in  megahertz).  The  decorrela¬ 
tions  due  to  the  sea  motion  for  the  frequencies  of  interest  are  almost  the  same  as  those 
given  in  Table  6  for  a  2-degree  beamwidih.  Consequently,  the  2-degree  beamwidth  de- 
correlations  should  be  used  for  any  beamwidths  smaller  than  2  degrees. 


Table  6 

Decorrelation  Time  for  A  Forward-Scanning  Radar 


Frequenev 
(.MHz)  ■ 

/  (sec)  for  Various  Beannvidtlis  i 

0.5 

1.0 

!  1.5 

i 

2.0  J 

140 

13.869 

3.467 

,  1.541 

0.866  1 

220 

8.826 

2.206 

I  0.980 

0.551  1 

440 

4.413 

1.103 

i  0.490 

0.275  ' 

900 

2.157 

0.539 

!  0.239 

0.134  . 

1300 

1,493 

0.373 

1  0.166 

0.093  j 

2900 

0,669 

0.167 

'  0.074 

0.041  1 

5290 

0.369 

0.092 

,  0.041 

0.023  1 

8500 

0.228 

0.057 

:  0.025 

o 

o 

Since  the  decorrelation  times  arc  still  quite  largo,  the  integration  gain  associated 
with  sea  clutter  will  be  less  than  tlie  integration  gain  associated  with  the  noise.  When 
the  sea  clutter  samples  are  correlated,  the  integrated  sigiial-to-cluttor  plus  noise  ratio 
should  l:ie  calculated  as  follows:  First,  the  signal  (clutter  +  noise)  ratio  lor  one  pulse  is 
given  by  >  ■/  .  Ni  ,  and  the  ratio  lor  r  pulses  is  given  apiiroxiniatcly  by 


;  1-  .  111 

where  /  ■  V‘  is  llie  integialion  gain  normally  associated  w  ith  pulses  and  C  is  the  effec¬ 
tive  number  of  independeiu  pulses  of  sea  duller  integrated: 
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Waves  Generated  by  a  Local  Wind  at  Sea 
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The  i>rf  is  the  pulse  repetition  rate  of  the  radar. 

The  question  remains  as  to  how  much  one  must  squint  the  beam  before  all  the  pulses 
will  decorrelate.  Since  tlie  unambiguous  range  of  a  forward-scanning  radar  requires  a 
pr(  of  about  65,  a  decorrelation  time  of  about  0.015  .second  is  necessary.  T  us,  using 
Table  8,  which  gives  the  decorrelation  time  at  a  4-degree  squint  angle,  one  can  see  that 


Table  8 

Decorrelation  Time  for  a  4-degree  Squint  Angle 


Frequency 

/'  (sec)  for  Various  Beamwidths 

(MHz) 

0.5 

- 

1.5' 

2.0' 

140 

0.446 

0.216  ! 

0.139 

0.101 

220 

0.284 

0.137  ] 

0.089 

0.064 

440 

0.142 

0.068 

0.044 

0.032 

i  900 

0.069 

0.033  1 

0.021 

0.015 

1300 

0.048 

•  0.023  ! 

1 

0.015 

0.011 

2900 

0.021 

0.014  ; 

0.006 

0.005 

5250 

0.011 

0.005  ' 

0.003 

0.002 

1  8500 

0.007 

0.003  1 

0.002 

0.001 

one  will  luive  couqjlete  independence  of  the  clutter  samples  if  one  uses  any  oi  the  follow¬ 
ing  antennas  and  frequency  combinations: 

2.0-degree  antenna  at  900  MHz 
1.5-degree  antenna  at  1300  MHz 
1.0-degree  antenna  at  2900  MHz 
0.5-degiee  antenna  at  5250  MHz  and  6500MHz. 

Thus,  if  one  has  any  of  the  comliinations  given  above,  one  will  have  complete  independ¬ 
ence  of  clutter  outside  the  4-degree  sector  and  partial  dependence  ot  clutter  inside  the 
4 -degree  sector. 

If  tlie  power  necessary  to  detect  a  given  size  target  along  the  scanned  area  were  cal¬ 
culated,  an  exorbitant  amount  of  power  would  be  required  in  Die  center  section  due  to  the 
clutter  correlation.  Thus,  to  reduce  tlie  required  power,  it  seems  desirable  to  slow  down 
the  scan  rate  in  the  center  section  so  as  to  obtain  more  independent  clutter  sampie.s.  If 
one  is  completely  clutter  limited  (no  tliermal  noise),  an  estimate  of  Die  slow-down  rate 
can  be  outaiiied  liy  dividing  the  desired  decorrelation  lime  (0.015  second)  by  the  decorre¬ 
lation  times  of  the  sea  clutter.  For  instance,  at  1300  MHz,  Table  6  gives  a  decorrelation 
time  of  0.093  second.  Thus,  the  slow-down  rate  should  be 
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^  =  0.16  ,ii  1-100 

MIb 

(6) 

0  003 

^  0.  t:  Ul  2900 

0.  oil 

Mil/  , 

(7) 

Since  one  is  not  really  clutter  limited,  the  optimal  slow-down  rates  will  be  larger 
than  those  given  in  Eqs.  (6)  and  (7).  In  fact,  a  slow-down  rate  of  0.5  is  probably  near  the 
optimum. 


TARGET  FLUCTUATIONS  (U) 

Point  Scatterers  as  a  Target  Model  (U) 

Most  targets  can  be  considered  to  consist  of  many  individual  point  scatterers.  If 
one  would  look  at  the  target  from  different  aspect  angles,  the  target's  cross  section 
would  appear  to  fluctuate  because  of  the  cancellations  and  reinforcements  between  the 
different  scatterers. 

To  determine  what  the  target  model  should  be,  an  approach  similar  to  one  used  in 
calculating  the  sea  clutter  decorrelation  time  was  adopted.  First,  from  the  geometry 
(Fig.  2)  one  can  write 

-  /  -  I  Cos  -t  . 
r  =  /  s  I  n  :  , 

^  y-  -  -  .  (8) 

<1!  •  ■’  Jl'i-  =  !  tj  -  !  -,0  -  -  (  X  -  I  ■  .  (9) 


Fi2,  Z  -  Radar  pcoinetry  for  t.irpet  a.i - 

tions,  with  the  largft  oi  longtl-,  !  shown  at 
one  ext  re  me  of  the  beam  width 
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Subtracting  Eq.  (8)  from  Eq.  (9)  one  obtains 

211-'. K  t  '  R-  --  <  ’.y-  2r'  .r  -  '  r'  . 

Consequently,  the  one-way  range  difference  is  approximately 

.  ./?  '  I  V  y  T  z  -  -x)  a  - 


and  the  tw'o-way  phase  change  is 


where  is  the  radar  wavelength.  Usually,  the  number  of  independent  samples  as  the 
beam  sweeps  past  the  target  is  taken  to  be  equal  to  the  number  of  phase  shifts  across 
the  beam.  At  one  extreme  of  the  beam,  ,  -  nt  2  and 


(10) 


At  the  other  extreme,  ■  -  n  2  and 


Subtracting  Eq.  (11)  from  Eq.  (10),  dividing  by  :  and  adding  unity  yields  the  number  of 
independent  samples: 


I  n  'll- 1' .  -  - ~  (  I 

The  returning  signal  from  \  equal  scatterers  is 


(12) 


or,  equivalently, 


-  s  i  n  ■ .  :  'y  cos  ; ,  -  cos 

and  the  cross  section  of  the  target  ,  is  proportional  to 


While  the  target  passes  through  the  beam,  the  target  will  appear  to  fluctuate  because 
of  the  changing  phases.  Typical  fluctuations  are  shown  in  Fig.  3  for  two  target  orienta¬ 
tions.  For  the  ;  =  0.0-radian  case,  Eq.  (12)  indicates  the  presence  of  only  2.7  independ¬ 
ent  samples;  and  for  the  ^  =  1.0  case,  Eq.  (12)  indicates  1.9  independent  samples.  If  one 
would  use  the  common  method  of  block  correlation  with  this  number  of  independent  sam¬ 
ples,  there  would  exist  a  probability  of  obtaining  a  very  low  strength  signal  requiring  a 
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F;c.  3  -  Ri;t  ii'ncd  sigr.al  irom  a  I’iutluatir.g  target 
{!  =  100  It;  =  1.0  It;  \  =  B;  ••  =  0.5  degree) 


very  high  transmitted  power  to  detect  it.  However,  if  one  took  a  very  large  number  of 
samples  (large  in  relation  to  the  number  of  independent  samples),  one  would  essentially 
obtain  the  average  returning  signal  over  the  given  aspect  angle;  and  it  appears  from 
Fig.  3  that  because  of  the  correlation  between  samples  one  will  not  obtain  the  small 
cross  sections  possible  when  using  block  correlation  in  connection  with  Eq.  (12).  To  test 
this  conjecture,  the  following  experiment  was  conducted;  For  each  of  100  cases,  eight 
scatterers  were  uniformly  distributed  over  a  length  of  100  feet  oriented  at  an  angle  of 
:  =  0.0  radians.  For  each  distribution,  tiie  average  returning  signal  s  was  found. 

Then,  the  average  and  standard  deviation  of  the  100  '.s,  were  calculated: 

-  J.  3 

and  (14) 

1  \  ,if  >  .  !'  -  .  0  'i3  . 

Since  what  one  is  really  trying  to  do  with  block  correlation  and  independent  samples  is  to 
model  the  actual  situation,  wliat  one  wants  to  find  is  the  number  of  independent  samples 
necessary  to  approximate  the  actual  distribution  —  or  equivalently  in  this  case,  the  stand¬ 
ard  deviation.  To  accomplish  this,  in  the  beginning  100  random  samples  were  taken 
from  the  100  configurations,  and  100  averages  v  each  consisting  of  only  one  point  were 
found.  Then,  the  standard  deviation  of  the  100  \,  •  was  calculated: 

■  V \  ‘  ;  I  j  , 
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Since  1.42  is  greater  than  the  true  standard  deviation,  more  than  one  independent 
sample  exists.  To  find  the  exact  number,  the  above  procedure  was  repeated.  However, 
this  time,  each  A,  consisted  of  two,  three,  and  four  independent  samples.  The  standard 
deviations  of  these  cases  are  as  follows: 


A,  consists  of  two  samples. 

A,)  : 

'  -  -  0.  'JH  . 

(15a) 

A,  consists  of  three  samples. 

;wir>  A,  I ) '  -  0.  :r  . 

(15b) 

A  consists  of  four  samples. 

,  V  a  r  ‘  \  )  . 

'  -  -  061. 

(15c) 

Thus,  comparing  Eqs.  (15)  to  Eq.  (14)  it  appears  that  there  are  approximately  four 
independent  samples  present  compared  with  the  2.7  indicated  by  Eq.  (12).  Repeating  this 
procedure  several  other  times  always  yielded  the  same  result  —  the  number  of  independ¬ 
ent  samples  required  to  equate  the  standard  deviations  was  always  slightly  larger  than 
the  number  given  by  Eq.  (12). 

In  the  rest  of  this  section,  the  number  of  independent  samples  will  be  defined  to  be 
that  number  for  which  the  standard  deviations  of  A'  and  >  are  equal.  Before  proceeding 
to  calculate  the  number  of  independent  samples  as  a  function  of  aspect  angle,  the  problem 
will  be  reformulated  to  include  ship  roll. 


Reformulation  to  Include 
Ship  Roll  (U) 

The  coordinate  systems  for  the  radar 
platform  and  for  the  target  are  shown  in 
Fig.  4.  The  following  notation  and  defini¬ 
tions  will  be  used; 

.■■i  axis  which  is  the  line  on  the  earth's 
surface  perpendicular  to  the  radar 
track, 

.  axis  which  is  the  line  on  the  earth's 
surface  parallel  tothe  radar  track, 

,  axis  which  completes  the  right- 
handed  coordinate  system, 

/.  iixis  which  is  the  line  from  the 
radar  to  the  target, 

axis  which  is  perpendicular  to  . 
and  lying  in  the  .  , /  .  plane. 


RADAR 

TRACK 


^  lu*  4  -  a !'.d  ldi‘u^‘1 

coorciir.  ;ii'  ivsK  ins 


axis  which  completes  the  right-handed  coordinate  system, 


squint  angle  of  the  radar,  i.e.,  the  angle  from  the  radar  to  the  target  center 
measured  with  respect  tothe  .  axis  in  the  ■  plane, 

angle  the  center  line  of  the  target  makes  with  the  -  ,  axis, 

roll  angle  of  Hie  target, 

grazing  angle. 
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Let  us  initially  assume  the  target  is  aligned  along  the  axis  =  0)  and  the  prob¬ 
lem  is  to  find  the  coordinates  of  a  point  in  the  y  system  whose  coordinates  in  the  .  sys¬ 
tem  are  u, ,  ,  and  w,.  First,  the  new  coordinates  due  to  a  roll  angle  are  given  by 


.  I 


1 

I 


SI 


i  *  ° 

.  0  rns 

I 

jo  sin 


0 


-s  I  n 


c  (1  s 


(16) 


To  line  up  the  j-,  and  y  .  axes  and  to  compensate  for  the  fact  that  the  target  is  not 
aligned  along  the  axis,  it  is  necessary  to  rotate  about  the  j ,  axis  an  angle  of  (  -  ); 


,  /  I  '  ‘  I  1  *>  I 

I  I  I  ■ 

Finally,  one  needs  to  rotate  through  the  grazing  angle  : 

I  0  -  ^  n:  j  '  ^  ‘ 

j  0  1  0  I.'  - 

I  >  i  11  ■  n  i-  r  1  M  I  <  ,  , 


(17) 


(18) 


Points  with  the  same  coordinate  are  at  essentially  the  same  range.  Hence,  the 
relative  phase  change  with  respect  to  the  center  of  the  target  is  given  by 

■...-5--./,  (19) 

Substituting  Eqs.  (16),  (17),  and  (18)  into  Eq.  (19)  yields 


(20) 


L-  J 

Next,  Eqs,  (13)  and  (20)  were  used  in  a  Monte  Carlo  investigation  to  determine  the 
number  of  independent  samples  for  side-looking  and  forward-scanning  systems  as  a 
function  of  the  orientation  of  the  target.  In  lieu  of  the  absence  of  a  standardized  sliip 
model,  the  following  ship  parameters,  which  do  not  seem  unreasonable,  were  used; 


Length  of  target,  150  feet, 
Width  of  target,  20  feet. 
Height  of  target,  13  feet. 
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Roll  angle  sin  ?  J- .■  7,  -.  i, 

n.uv>  0-2  radian, 

7p ,  10  seconds. 

Also,  the  following  radar  parameters  were  used: 

Beamwidth,  0.5  degree. 

Grazing  angle,  3.5  degrees. 

Integration  time,  1.5  sec  for  a  for^’ard-scanning  system, 

3.2  sec  for  a  side-looking  system. 

For  each  of  the  100  configurations  used,  14  scatterers  were  uniformly  distributed  in 
the  rectangular  target  volume,  and  the  initial  phase  of  the  roll  .  was  uniformly  distrib¬ 
uted  on  0  to  J  .  The  number  of  independent  samples  tor  the  side-looking  radar  is 
shown  in  Fig.  5  and  the  number  of  independent  samples  for  the  forward-scanning  radar 
is  shown  in  Fig.  6.  For  the  side-looking  radar,  the  number  of  independent  samples  in¬ 
creases  from  one  to  24  as  the  target  orientation  changes  from  bow  to  broadside;  and  for 
the  for\^■ard-scanning  radar  the  variation  is  one  to  11  from  bow  to  broadside.  To  see 
how  much  these  fluctuations  cost  w-ith  respect  to  complete  independence  between  pulses, 
one  can  consult  Table  9,  which  was  derived  from  Swerling's  report  (4)  about  log-normal 
fluctuating  targets. 


Table  9 

Required  Increase  in  the  \  Ratio  with  Respect  to  Complete  Independence' 


Probability 
of  Detection, 

/■' , 

s  N  Ratio  Increase  for  Decreasing  Numbers  I 

of  Independent  Samples  (db)  | 

100 

10 

4 

^  .  i 

.99 

.9 

0  !  2 

0  1.25 

4.75 

3 

7.25  11.75 

4.5  7.5 

_ ; _ 

"Thu<(;  for  a  mear. -to -media r,  raf.o  of  Z.  a?  one  (.an  sc’(; 

from  G'lari;  (Ugiii'.i  (S)  cl'.ese  nonibor?  c.xe.  be  ,.se(i  lor  a  R.'yieijh  fine;  ;aling 
t.-i-.-gei. 


If  one  wanted  id  maintain  a  probability  of  detection  of  .99  for  all  target  orientations, 
It  would  be  necessary  to  increase  the  5  A  ratio  by  either  7.25  or  11.75  db,  depending  on 
the  operating  frequency.  This  seems  like  an  exorbitant  waste  of  power  which  is  not  jus¬ 
tified  because  of  the  following  discussion. 

Even  though  the  probability  of  detection  is  greatly  deteriorated  at  bow  aspects,  it  is 
very  little  deteriorated  at  broadside.  A  natural  question  to  ask  is  what  is  the  overall 
probability  of  detection  assuming  that  ships  will  be  uniformly  distributed  in  aspect.  By 
using  Swerling's  curve  (Fig.  7)  and  Fig.  5,  a  side-looking  radar  operating  at  =  0.333 
feet  has  the  following  minimum  probabilities  of  detection: 
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0  0 

<  < 

0.  1  . 

mi  n  ) 

-0.64 

0.  1 

<  •  < 

0.  t  . 

min  *. /* ;  .' 

=  0.76 

0.  i 

—  -> 

mi  n  * j  ' 

--  0  00 

Hence,  the  overall  probability  of  detection  Is 

r.  =  !o.  1  i  (I- ()  t)  i  0.  .<'■  0.  re,)  •  '  '2-0  dio.O''  -  . 

The  calculation  was  performed  for  various  integrated  \  ratios,  frequencies,  and  the 
two  radar  configurations.  The  results  are  shown  in  Table  10. 


Table  10 

Probability  of  Detection  for  Fluctuating  Targets' 


Integrated 
ft  V  Ratio 
(db) 

Side-Looking  Radar 

Forward-Scanning 

Radar 

=  1  ft 

=  0.33  ft 

■  -  1  ft 

=  0.33  ft 

}  =  0.2  ft 

16 

.72 

.86 

.62 

.73 

1  .76 

17 

.81 

.91 

.72 

.82 

.85 

18 

.88 

.96 

.80 

.89 

!  .91 

19 

.92 

.98 

.86 

.93 

I  .95 

20 

.95 

.99 

.91 

.96 

1  .97 

’  Ar,  iritk'iir.iU'cl  >  V  ri'.lio  of  lofi'')  <. o r '.'C' sponds  to  probability  of  false  .ilariri  of  10*''  ar.d  a 
prai);i  hil  Uv  of  detoi't  ioo  of  .9-'’  for  I  nr  -,is(  \'.in  ,■•!!  Ibc  r,,  irt-  i  nfi  ciu  I'.dei'it . 


Using  linear  interpolation  on  Table  10,  the  following  r  \  ratios  are  required  for  a 
probability  of  detection  of  .9  for  fluctuating  targets  and  a  probabilily  of  false  alarm  of 


10'“' .  For  a  side-looking  radar 

\  in.  I.  .11.  (.,r  •  -  n.  i.l  ii  ,  (21a) 

>  \  I  t  ,ii.  r...  ,..i  fi  .  (21b) 

For  a  forward-scanning  radar 

\  i:  K  .11.  i..i  (1  2  fi  ,  (21c) 

>■  \  lit.:  .11.  I  0  1 1  ,  (21d) 

'  \  |i).  'I  .11.  i...  ■  I .  n  ( I  .  (21e) 


Thus,  witli  the  integrated  ratios  shown  in  Eqs.  (21),  one  obtains  a  probability  of  de¬ 
tection  of  over  .9  for  fluctuating  targets  and  a  probability  of  detection  of  over  .99  for 
nonfluctuating  targets. 
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The  ■?  \  ratios  in  Eqs.  (21)  are  those  required  vhen  the  target  is  seen  at  a  grazing 
angle  of  3.5  degrees  for  3.2  sec.  If  the  grazing  angle  is  increased,  the  decorrelation 
times  due  to  rolling  will  change,  and  new  s  \  ratios  will  be  required.  If  the  grazing  an¬ 
gle  is  increased  to  51  degrees  with  a  corresponding  integration  time  of  0.44  sec,  the 
number  of  independent  samples  for  a  side-looking  radar  is  shown  in  Fig.  8.  Comparing 
Figs.  5  and  8,  one  sees  that  one  has  more  independent  samples  at  bow  aspect  but  less 
independent  samples  at  broadside  for  the  larger  grazing  angle.  Using  the  same  philos¬ 
ophy  as  before,  the  probabilities  of  detection  have  been  recalculated  for  the  larger  graz¬ 
ing  angle  and  are  shown  in  Table  11.  The  values  in  Table  11  were  calculated  using  Fig.  9 
instead  of  Fig.  7,  because  at  the  large  grazing  angles  (small  ranges)  one  is  integrating 
only  29  pulses  in  comparison  with  over  100  at  the  large  ranges. 


FTc.  S  -  I  Ml-  1-  of  ir.di  pn-dctu  sample;, 

for  .1  side- lookui,;  radar  (pra/ir.e  ar.jjir  =  51 
r«-i  ;  iiiUgratior.  liinc  =  O.t-i  ifC'C) 


Table  11 

Probability  of  Detection  for  Fluctuating  Targets  (Grazing  angle  =  51 
degrees;  Integration  Time  -=  0.4  sec;  Side-Looking  Radar) 


Integrated 
'  V  Ratio 
(db) 

I’. 

Integrated 
\  Ratio 
(db) 

■  =  1  ft  !  =  0.33  ft 

■■  =  1  ft  1  =  0.33  ft 

16 

17 

.80  ;  .88 

.86  I  .92 

18 

19 

.91  '  .96 

i 

.94  1  .98 
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Fig.  ^  -  Radar  deti;ction  probabii- 
itv  lor  log  -  iic  rcTia  il  V  distrib.J.ed 
tarucds  (  =  2,  ^  -  20;  probabililv 

of  a  false  alarm  =  10  } 


Again  using  linear  interpolation  on  Table  11,  the  >  A  ratios  required  for  a  probabil¬ 
ity  of  detection  of  .9  and  a  probability  of  false  alarm  of  lO"' "  with  a  grazing  angle  of  51 
degrees  and  an  integration  time  of  0.44  sec  for  a  side-looking  radar  are 


s  \  ! :.  !i  -11.  f.M  ■  : . "  f  I  .  (22b) 

As  one  can  see,  these  values  are  slightly  smaller  than  those  given  in  Eqs.  (21). 

Since  the  values  in  Eqs.  (22)  were  calculated  using  Fig.  9,  which  assumes  20  pulses  in¬ 
tegrated  whereas  really  29  pul.ses  are  integrated,  the  required  >  \  ratios  will  be  slightly 
larger  than  those  given  in  Eqs.  (22).  Consequently,  the  desired  >  v  ratio  can  be  consid¬ 
ered  a  constant  over  all  grazing  angles. 

The  required  v  A  ratios  have  br-rn  calculated  for  several  otl'er  integration  times. 
The  results  of  the  various  calculations  are  presented  in  Figs.  10,  11,  and  12. 


IONOSPHERIC  EFFECTS  (U) 

This  section  discusses  how  the  ionosphere  affects  detection  of  targets  at  sea.  The 
two  major  problems  caused  Ijy  passage  of  the  signal  through  the  ionosphere  are  the  Far¬ 
aday  rotation  of  the  signal  and  the  random  phase  shift  associated  with  the  received  sig¬ 
nal."  The  approach  consist.^  of  postulating  a  very  simplified  model  of  llicse  effects  and 
then  calculating  the  degradation  caused  by  these  effects. 
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Fig.  12  -  Ef;%'ct  ot  tl-.e  integration  lime  on  the  rc- 
quiVvcl  ratio  lor  a  forward-ecanr.ing  radar  (ire- 

qacni  y  =  1300  Me;  grazing  angle  =  2  degrees) 


Faraday  Rotation  (U) 

Faraday  rotation  is  a  complicated  function  of  geographical  position,  time  of  day, 
time  of  year,  orientation  of  the  radar  beam  with  respoct  to  the  earth’s  magnetic  field, 
and  many  other  factors.  Therefore,  instead  of  forming  an  extremely  complicated  model 
for  Faraday  rotation,  the  postulated  model  will  only  reflect  some  ot  the  worst  conditions 
to  be  encountered:  upper  midlatitudes,  noon  wintertime,  20-degree  depression  angle, 
and  orientation  parallel  to  the  earth's  magnetic  field.  For  these  conditions,  by  assuming 
a  Chapman  distribution  as  a  model  for  the  ionosphere  electron  density  distribution,  the 
Faraday  rotation  can  be  approximated  (6)  by  the  following  equation  and  table; 

=  •  lo  '  f-  .  (23) 


Altitude  II 
(naut  mi) 

'■)/ 

Altitude  II 
(naut  mi) 

- — - J 

150 

1.0 

109 

300 

0.945 

322 

200 

0.97 

198 

400 

0.92 

377 

250 

0.956 

272 

600 

0.91 

413 
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where  t-  is  the  one-way  Faraday  rotation  in  radians,  /  is  the  frequency  in  hertz,  Cn  is 
related  to  the  strength  of  the  magnetic  field  and  is  a  function  of  the  altitude  //,  and  Cu  is 
a  slab  thickness  parameter  and  is  a  function  of  the  altitude  II.  Using  Eq.  (23)  and  the 
table,  the  one-way  Faraday  rotations  were  calculated  and  are  shown  in  Table  12. 


Table  12 

One-Way  Faraday  Rotation  of  the  Signals  at  Various 
Frequencies  With  the  Satellite  at  Various  Heights 


Altitude  II 
(naut  mi) 

Faraday  Rotation  ■'  (radians) 

140 

MHz 

220 

MHz 

440 

MHz 

900 

MHz 

1300 

MHz 

2950 

MHz 

5250 

MHz 

8500 

MHz 

150 

12.3 

5.0 

1.25 

0.30 

0.14 

0.03 

0.01 

0.00 

200 

21.7 

8.78 

2.20 

0.53 

0.25 

0.05 

0.02 

0.01 

250 

29.4 

11.9 

2.98 

0.71 

0.34 

0.07 

0.02 

0.01 

300 

34.4 

13.9 

3.49 

0.83 

0.40 

0.08 

0.02 

0.01 

400 

39.2 

15.9 

3,96 

0.95 

0.46 

0.09 

0.03 

0.01 

600 

42.5 

17.2 

4.30 

1.03 

0.49 

0.10 

0.03 

0.01 

_ 

To  calculate  the  degradation  in  the  signal-to-clutter  ratio  and  .s  \  ratio,  it  is  neces¬ 
sary  to  calculate  the  received  signal  after  it  has  passed  through  the  ionosphere,  been  re¬ 
flected  at  the  sea  surface,  and  again  passed  through  the  ionosphere.  First  of  all,  let  us 
assume  that  the  transmitted  electric  field  is 


E-  =  ,E:<  .  1  -■  -s  ' 

-  E,'//-yV‘l  . 

where  y  is  the  cross-polarization  gain  of  the  antenna  and  n  and  ''  are  unit  vectors  in  the 
horizontal  and  vertical  directions.  Then,  the  electric  field  reflected  by  the  sea  is  given 
by 

=  6,  \  'll  ’,<■(»>  :  r  't  -  </  sin  '  •  ^  ill  ^  ‘  \  (-  'si  n  <  ■  -  y  ct's  ■  +’!'(  ,  (24) 

where  is  the  sea  clutter  cross  section  for  horizontal  polarization,  '■  is  the  sea  clut¬ 
ter  cross  section  for  vertical  polarization,  is  the  one-way  Faraday  rotation,  r  is  the 
phase  shift  between  the  horizontal  and  vertical  components  that  is  caused  by  the  reflec¬ 
tion,  and  '  is  the  incident  angle  between  the  horizontal  component  of  6'.  and  the  sea  when 
there  is  no  Faraday  rotation.  In  reality,  -  is  a  random  variable  caused  by  the  roughness 
of  the  sea.  Since  its  distribution  is  unknown  and  since  it  is  not  of  prime  importance  in 
this  discussion,  we  will  assume  -  o  .  Then  Eq.  (24)  reduces  to 

-  El  ‘  I  .IS  -  ‘1  sjp  ill  ‘  f  '  \  j  ■■  s  1 11  ■  ■/  ■  iis  .  1  I 

The  power  received  by  an  antenna  is  proportional  to  the  square  of  the  absolute  value 
of  the  received  electric  field.  That  is 
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- <'  *  *  ^  -  s  i  n  ^  •  +  ij 

-  n 


^  t  f)  3  “  :  ■  I 
V‘  cos‘:)| 


represents  the  backscatter  from  sea  clutter.  On  the  basis  of  work  by  Macdonald  (7,8)  a 
ratio  of  10  db  was  selected  as  a  good  average.  Also,  it  is  fairly  easy  to  build  a 
linearly  polarized  antenna  whose  cross-polarization  response  is  down  at  least  20  db 
(;,•  =  0.1).  Finally,  while  .  varies  over  a  wide  range  of  values,  ,  =  n  was  chosen,  since 
this  value  yields  the  tvorst  losses.  For  these  values  Eq.  (25)  gives  for  the  backscatter 
power  from  sea  clutter 


A  similar  approach  can  be  used  to  find  the  returning  signal  power;  the  only  differ¬ 
ence  is  ihat  the  ratio  ■  ,  of  the  target  cross  sections  for  vertical  and  horizontal 
polarizations  (again  based  on  Refs.  7  and  8)  is  -4  db.  Thus,  the  returned  signal  power  is 


The  bracketed  term  in  Eq.  (27)  represents  the  degradation  of  the  >  ^  ratio,  and  the 
ratio  of  Eq.  (27)  and  Eq.  (26)  represents  the  degradation  of  the  signal-to-clutter  ratio. 
These  losses  are  presented  in  Talile  13. 

Table  13 

Degradation  of  the  >  Ratio  and  Degradation  in  the 
Signal-to-Clutler  Ratio  Caused  by  Faraday  Rotation 


Altitude  /.' 
(naut  mi) 


Degradation  Factors  at  Various  Frequencies 

;20  •  440  I  900  I  1300  I  2950  1  5250  I 

;H2  '  MHz  i  MHz  i  MHz  I  MHz  '  MHz  i 


8500  i 
MHz  I 


Signai-to-Noise  Power  Degradation  Factor  (Factor  in  Eq.  (27)) 


'  0.0001 
'  0.0001 
!  0.0001 
i'  0.0001 
]  0.0001 
!  0.0001 


0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 


0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 


0.743  I 
0.340  I 
0,094  ' 

0.0120  ! 
0.0001  j 
0.0001  I 


Signal-to-Clutter  Degradation  Factor  (Ratio  of  Eq.  (27)  and  Eq.  (26)) 


I  0.0001 
,  0.0001 
!  0.0001 
'  0.0001 
i  0.0001 
;  0.0001 


0.0001 

0.0001 

0.0001 

0,0001 

0.0001 

0.0001 


0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 


1.0 

1.0 

0.1632 

0.0077 

0.0001 

0.0001 


1.0 

i  1.0 

1  1.0 

1.0 

1  1.0 

1.0 

i  1.0 

1  1.0 

1.0 

1  1.0 

!  1.0 

1.0 

1  1.0 

1  1.0 

1.0 

1.0 

!  1.0 
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In  deriving  Table  13  the  following  method  was  used;  For  a  given  frequency  and  alti¬ 
tude,  the  losses  recorded  in  Table  13  correspond  to  the  maximum  losses  associated  with 
any  rotation  between  zero  and  tlie  maximum  Faraday  rotation  that  is  given  in  Table  12. 
This  is  why  the  maximum  loss  of  0.0001  is  so  frequent.  Also,  lor  small  rotations  the 
ratio  of  Eqs.  (27)  and  (26)  is  greater  than  1,  since  no  depolarization  has  been  assumed; 
i.e.,  -■  .  Consequently,  a  conservative  approach  was  taken  by  replacing  these  gains 

by  I's. 

It  seems  likely  that  one  could  reduce  the  losses  for  the  low  frequencies  il  one  would 
use  circular  polarization.  Tliis  is  equivalent  to  assuming  that  - .  Then,  Eqs. 

(26)  and  (27)  reduce  to 

fi  /  _ 

'  6-  ,  (^1  ■  N  .  I  :  . 

and 


Thus,  the  \  power  degradation  factor  is  0.068  and  the  signal-to-clulter  power  degra¬ 
dation  factor  is  0.0297  for  circular  polarization,  these  factors  being  independent  of  the 
Faraday  rotation. 


Random  Phase  Siiift  (C) 

Another  effect  of  tlte  ionosphere  is  the  random  pliase  shift  introduced  by  variations 
in  the  iou  density.  This  phase  shift,  whose  mean  is  directly  proportional  to  time,  limits 
the  coherent  integration  time.  Consequently,  the  question  arises  as  to  how  many  pulses 
\  one  sliould  integrate  coherently  if  one  receives  v  pulses  in  time  /  .  Before  calculating 
the  optimal  coherent  integration  gain  when  the  mean  phase  of  the  received  signal  drifts 
linearly  with  time,  the  total  integration  gain  will  be  calculated.  First,  the  received  cn- 
erg>’  for  v  pulses  is  given  liy  the  correlation 


(28) 


where  is  the  time  between  pulses  =  i  |.ri  ,  is  tlie  pulse  width,  and 
dom  phase  shift,  which  will  be  approximated  by  its  mean;  that  is  ■  ■  ■ 
the  phase  shift  per  unit  time  and  therefore  is  a  function  of  the  frequency 
Eq.  (28)  yields 


is  the  lan- 
whorc  is 
Integrating 


(29) 


which  for  the  's  of  interest  can  be  approximated  by 


(3U, 
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Expanding  Eq.  (30), 


i'  = 


(31) 


and  substituting  the  small  angle  approximations 


1 


into  Eq.  (31),  one  obtains 


V-  I 

....  .-,r. 


(32) 


It  can  be  shown  that  the  summation  in  Eq.  (32)  can  be  rewritten  as 


■  - 


...  b  -  AK  r] 


>  I  „  ■  ;  7 . 


which  if  one  uses  the  small  angle  approximation  reduces  to 


i  I  '  j  •  I 

^  T[T  =  - 


-■  r. 


or 


=  T 


,  s  I  11  (  A  -  4)  •'  I  .. 

^  .  .1)  -Atv-'- 


If  V  is  large,  this  can  simply  be  approximated  by 


■  '  s,„  'l>\ 

^  ■  7  I'  —!  ■ 


(33) 


where 't  is  the  phase  shift  encountered  in  time  /  .  •  \  -  i  -  >  : 

I  -  -t''  • 

The  term  in  parentheses  in  Eq.  (33)  represents  the  coherent  integration  gain  for  \  pulses 


\  ...i  'I- 
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Thus,  of  the  w  pulses  one  receives,  N  consecutive  pulses  are  coherently  integrated;  and 
then,  the  results  of  these  coherent  Integrations  are  combined  by  performing  a  noncoher¬ 
ent  integration  on  A  =  v  v  pulses.  The  noncoherent  integration  gain  is  given  by 


S(A)  c 

1 

fo  r 

A  =  1  . 

(34a) 

AlA)  c 

1  .  0 1  0  A  •  •  ’  “ 

f  o  r 

1  <  A  <  t  . 

(34b) 

S(A)  = 

CO 

f  o  r 

■1  1  A  ■-  20  . 

(34c) 

■SIA'i  - 

1 .  6t:>  k-' 

for 

20  1  A  <■  100  . 

(34d) 

A'AI  c 

"94  K  -  ■ 

f  o  r 

A  100  . 

(34e) 

which  are  curves  that  Fox  (9)  has  used  to  approximate  the  noncoherent  integration  gain 
given  in  Skolnik  (10).  Therefore,  the  gain  associated  with  integrating  t'  pulses,  \  pulses 
coherently,  is 


where 


A-  . 


(35) 


The  problem  of  finding  the  optimal  integration  time,  or  equivalently  V,  is  reduced  to 
finding  the  value  of  v  that  maximizes  the  gain  given  in  Eq.  (35).  Setting  the  derivative  of 
Eq.  (35)  to  zero  yields  a  very  complicated  equation  for  the  optimal  value  of  \.  However, 
if  one  makes  the  approximation 


4>  -r  :r.  \  . 

the  equation  for  the  optimal  .\  is  quite  simple.  First,  <:  is  then  of  the  form 


Thus, 


Solving  for  \ , 


C  '■ 


s.n  .I'm/''' 


•  ? .  A 


!(■  I  i>  s  T  1 .  \  sin  i ^ 

(TV  ®  ^  \A-r 


;  7 ,  \  cm  ;  r_  S  ■  H  . 


(36) 


Consequently,  it  is  the  phase  shift  •  7  , \  which  limits  the  coherent  integration  time. 
The  solutions  of  Eq.  (36)  are 


^  0.41  . 

if  R  ^  0.014. 

(37a) 

4'  0.  80  . 

if  R  -  0. . 

(37b) 

O’  0-92. 

if  ■  0.  68.8  . 

(37c) 

4’  ■  I  .  06  . 

1  f  /£*  =  0.  .-,9.1  . 

(37d) 
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Thus,  the  limiting  phase  shift  is  seen  to  be  a  function  of  the  number  of  noncoherent 
pulses  that  are  integrated.  To  show  that  these  results  are  valid  in  spite  of  the  numerous 
approximations,  the  coherent  gain  equation,  Eq.  (29),  and  noncoherent  gain  equations, 
Eqs.  (34),  were  programed  on  the  CDC  3800  computer,  and  an  exhaustive  search  was 
used  to  find  the  optimal  value  of  v  and  ir  for  the  following  parametric  ranges: 

/  -  20,1  Mil,  . 


.  ^  I.  7122  . 
bO  |>  t  i  <210. 

I  2  1  -> . 

The  results  of  this  investigation  are  shown  in  Table  14.  One  can  infer  from  this 
data  that  when  one  is  integrating  a  short  time  (a  few  noncoherent  pulses)  the  limiting 
phase  shift  is  about  0.8  radian;  however,  when  the  integration  time  is  longer,  the  limiting 
phase  shift  is  increased  to  0.92  radian.  These  results  are  compatible  with  the  results  in 
Eqs.  (37).  The  loss  due  to  the  pliase  shift,  in  comparison  to  no  phase  shift  at  all,  is 
given  by  .  n  ,  which  is  0.9  or  -0.5  db  if  'i>  =  0.8  and  is  0.87  or  -0.6  db  if  =  0,9. 


Table  14 

Limiting  Phase  Shift 


prl 

\i  ■■  i-rr 

1/  =  1  |.i  f 

C  1-  I  1 

4' 

(radians) 

\' 

u _ 

■1. 

(radians) 

' 

i)> 

(radians) 

60 

11 

0.78 

■ 

11 

0.78 

12 

0.86 

80 

14 

0.77 

14 

0.77 

16 

0.88 

100 

18 

0.80 

18 

0.80 

20 

0.89 

21 

0.79 

21 

0.79 

24 

0.90 

24 

0.77 

24 

0.77 

28 

0,91 

160 

28 

0.80 

28 

0.80 

32 

0.91 

180 

31 

0.79 

31 

0.79 

36 

0.92 

200 

35 

0.80 

35 

0.80 

40 

0.92 

38 

0.79 

38 

0.79 

44 

0.92 

« 

0.79 

41 

0.79 

48 

0.92 

DATA  PROCESSING  (S) 

Optimal  Integration  Angle  (U) 

When  a  radar  scans  past  a  point  target,  the  amplitudes  of  the  returned  pulse  train 
are  proportional  to  the  square  of  the  one-way  antenna  voltage  gain  pattern.  Usually,  one 
considers  integration  over  the  beamwidth  ^  (which  is  detcrmint^  by  the  3-db  power 
points)  and  suffers  a  1.6-db  "scanning  loss,"  which  is  based  on  an  optimal  integration 
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angle  ol  0.84  ,  for  a  gaussian-shaped  beam.  However,  these  figures  are  based  on  the 

work  of  Blake  (ll),  who  assumed  a  noncoherent  gain  that  obeys  the  half-power  law.  To 
extend  the  work  to  other  laws,  a  slightly  different  approach  than  Blake's  will  be  used. 

Two  beam  shape  patterns  will  be  considered;  a  gaussiaii  beam  shape  and  a  -i  n  .  ■ 
beam  shape;  that  is, 


and 

-  ;  II  ■.  I  fl  (; 

I  >''!'■  - 

where  is  the  angle  measured  with  res[)ect  to  the  beam  center,  is  the  angle  of  the 
half-power  point,  and  the  gain  is  normalized  lo  unity  at  midbeam. 

To  find  tlie  optimal  integration  angle,  one  must  formalize  the  detection  problem  as  a 
test  of  a  binary  hypotliesis.  Therefore,  consider  the  problem  of  testing  for  a  known  sig¬ 
nal  in  white  gaussian  noise; 

II  /:■  .  ■  .  ■  .  \ .  .  u.  .  \ 

II.  :  I  t,  -t  1  .  - 

where  -'a  ■  i  is  the  number  of  pulses  integrated. 

The  null  liypothests.  .  is  the  statement  that  only  noise  is  present;  and  the  alterna¬ 
tive,  .'V,  ,  is  the  statement  that  signal  plus  iicise  is  present.  The  noise  samples  •' .  and 
are  independent  identically  distributed  gaussian  random  vai-iables  with  mean  zero 
and  \ariance  ■ .  Since  we  are  mostly  concerned  with  the  threshold  case  (small  >  A 
ratio),  the  square-law  detector  will  be  considered.  That  is,  the  lest  statistic 


z  ’ 


will  be  cmniijared  to  a  thrcsliold  / .  If  ilien  'V,  is  accepted  (signal  is  present'; 

and  if  v  /,  ihen  ll .  is  acceiitcd.  Since  is  equal  to  tlie  sum  of  independent  random 
variab'lcs,  is  essentially  gaussian  as  implied  by  the  central  limit  theurem.  Tiius, 
only  the  means  and  variances  of  under  boili  hypotheses  will  lie  imiioriant  in  the  de¬ 
tection  problem.  Consequently,  it  is  necessary  to  calculate  llieni.  In  the  forthcoming 
calculations  the  iollowing  well-known  information  about  Hie  moments  of  a  zero  mean 
gaussian  variable  will  be  used; 


SECRET 


32 


G,  V.  TRUNK 


SECRET 


-  Kvi/i')  ~  0  , 
h'l  j  ''  ^  =  3-*  . 


Then,  the  means  of  are 


IC-X'-lli  -  A 


r'  T  yf  I 


and 


-  A  I 

and  the  second  moments  are 


^  y  =  j  .  <v*,  *  2.-:  . 


I-:-/.'  II,)  A'  i-r  •  y-r  =  (■  -' 


and 


l-^/.f'll.l  -  .  y;|'  -  .)  »  H6”<^,)c•=  ‘  H-’  . 

The  variances  of  /;  are 

\.ir  >/,-  //.i  !■:</. -'ll'-  ■  Wv i»,l  •-  =4.' 

and 

vor:/2'//i  -  I  i:</.hii,': '  ■■  \  ' . 

Since  the  meaii  of  the  sum  of  random  variables  is  equal  to  the  sum  of  the  means  of 
the  random  variables,  and  since  tlie  variance  of  the  sum  of  independent  random  variables 
is  equal  to  the  sum  of  the  variances  of  these  random  variables,  one  obtains 


•  1- 


l■:■s^  ■■ll  >  ^  2  -  j'CN  :  I'  ■- 


(38) 


III, 


\  ’ 


\ 

\  ■ 


.1 


i 


(39) 


and 


. ,  ill  ' 


1  '  1 2  '  M 


(40) 


\ 

\  ’ 


Ah  I  T  1'  2\  I 


(41) 
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From  these  preliminaries,  one  can  proceed  to  the  heart  of  the  problem.  When  mak¬ 
ing  a  test  betu’een  two  gaussian  distributions,  r,-i and  (.'  i -  >  ,  the  probability  of 
detection  is  a  monotonic  function  of  the  expression  (as  proved  in  Appendix  A) 


I 


(42) 


where  'll  '  is  the  cumulative  distribution  for  a  normalized  gaussi?n  (i.e.,  oi  n  i ' ), 


ax  , 


and  consequently  "t"'! !  -  is  the  threshold  corresponding  to  a  probability  of  false 
alarm  of  l\  Using  the  appropriate  equivalences 

-  6''  ■ 


\ 


II,  ‘  . 


and 


:  V  -1  r  ^  S  \  /A, '  . 


ys.-  n  ■>' ' 

•  I 

substituting  Eqs.  (38)  through  (41)  into  Eq.  (42)  yields 


\  ^  A  .  K.  .  . 


<V  A  '  -  :  -  \:!N  ■  !  'I>'  I  -  l‘,  .) 


(43) 


where  one  receives  a  pulses  per  radian.  That  is, 

\  -  A.-. 

For  a  side -looking  radar  A  can  be  expressed  as 

,,  IllL'll 

’  i 

where  is  the  target  range  and  i  is  the  velocity  of  the  radar  platform. 

The  optimal  integration  angle  is  the  angle  that  maximizes  Eq.  (43).  It  is  quite 
obvious  that  the  optimal  ,  is  a  function  of  and  lienee  of  the  >  \  ratio.  However,  it  is 
surprising  to  note  that  tlie  optimal  ,  is  also  a  function  of  the  i>i  (  tlirough  the  variable  a  . 
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Equation  (43)  was  programed  on  the  CDC  3300  computer  with  <i>-'  {  \  j  =  6.232,*  „• 
varying  from  0.4  to  0.9,  and  K  varying  from  0.005  to  0.05.  The  results  for  a  i  sir  r'  j 
beam  shape  are  shown  in  Table  15. 


Table  15 

Optimal  Integration  Angle  for  a  isin  .r'  r  Beam  in  ,  Units 


— 

1  h 

■  (radians)  for  Various  \  Ratios 

4.95* 

3.01 

1.43 

0.09 

-1.07 

-2.10 

0.005 

1.145 

1.065 

1.005 

0.955 

0.915 

0.88 

0.010 

1.11 

1.03 

0.96 

0.910 

0.87 

0.84 

0.015 

1.08 

1.005 

0.93 

0.885 

0.84 

0.81 

0.020 

1.06 

0.98 

0.92 

0.86 

0.82 

0.78 

0.025 

1.05 

0.95 

0.90 

0.825 

0.80 

0.75 

0.030 

1.02 

0.93 

0.87 

0.81 

0.78 

0.75  i 

0.035 

1.015 

0.915 

0.84 

0.805 

0.77 

0.735 

0.040 

1.00 

0.915 

0.84 

0.80 

0.76 

0.72 

0.045 

0.99 

0.91 

0.81 

0.765 

0.72 

0.72 

0.050 

1.00 

0.90 

0.80 

0.75 

0.70 

o 

o 

■‘S  ,V  10  l,.,t  ■  1  -  •' 


The  question  naturally  arises:  If  one  integrates  over  the  optimal  angles  given  in 
Table  15  what  is  the  scanning  loss  associated  with  the  fact  that  the  returning  signals 
have  a  lower  strength  off  midbeam?  To  calculate  tliis  scanning  loss  it  is  necessary  to 
calculate  the  value  of  Eq.  (43)  for  a  rectangular  beamwidth.  Thus,  if  one  lets  n  ■  i , 
Eq.  (43)  becomes 


The  scanning  loss  for  a  particular  \  and  |.r(  will  be  defined  as  the  change  in  the  \ 
ratio  (i.e.,  the  change  in  )  required  to  equate  Eq.  (43)  to  Eq.  (44). 

To  illustrate  how  the  scanning  loss  is  calculated,  let  us  consider  the  case  where 
N  =  61.  Equations  (34)  give  an  integration  gain  of  14.47  db  for  an  integration  of  61 
pulses.  Consequently,  an  '  ratio  of  1.53  db  is  required  to  obtain  a  16-db  integrated 
.s  \  ratio.  Then,  Eq.  (43)  was  evaluated  tor  •  values  ranging  from  0.8  to  1.0  and  >  \ 
ratios  ranging  from  2.73  to  3.53.  The  results  are  presented  in  Table  16. 


‘'Till'  ill  risl.Qld  1 .)  ;■  ri' ipoiKi lo  /’,  ,  =  lO"' 
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Table  16 

Evaluation  of  Eq.  (43)  lor  v  =  61 


2.93 

3.43 

3.53 

0.310 

0.323 

,0.336 

0.349 

0.362 

0.375 

0.388 

0.401 

0.82 

0.290 

0.302 

0.315 

0.328 

0.340 

0.353 

0.366 

0.392 

0.84 

0.282 

0.294 

0.307 

0.319 

0.332 

0.345 

0.358 

0.384 

0.86 

0.274 

0.286 

0.299 

0.311 

0.324 

0.336 

0.349 

0.375 

0.88 

0.266 

0.278 

0.291 

0.303 

0.316 

0.328 

0.341 

0.353 

0.366 

0.90 

0.259 

0.271 

0.283 

0.295 

0.307 

0.320 

0.332 

0.345 

0.357 

0.92 

0.251 

0.262 

0.275 

0.287 

0.299 

0.311 

0.224 

0.336 

0.349 

0.94 

0.243 

0.255 

0.267 

0.279 

0.291 

0.303 

0.315 

0.328 

0.340 

0.96 

0.235 

0.247 

0.259 

0.271 

0.283 

0.295 

0.307 

0.319 

0.331 

0.98 

0.227 

0.239 

0.251 

0.263 

0.275 

0.286 

0.323 

0.220 

0.231 

0.243 

0.255 

0.266 

0.278 

■ 

0.314 

‘■.V  S  =  10  I  1 


Since  the  value  of  Eq.  (44)  is  0.313  for  \  =  61,  the  required  s  v  ratios  (which  can  be 
obtained  by  i  sing  linear  interpolation  on  Table  16)  are  shown  in  Table  17. 


Table  17 

Required  s  v  Ratio  for  Eq.  (43)  to  Equal  0.318  lor  =  51 


(radians) 

s  Ratio 

_ 

Scanning  Loss  (db) 

1  K 

0.80 

2.89 

1.36 

0.027 

0.82 

2.95 

1.42 

0.027 

0.84 

3.02 

1.49 

0.028 

0.86 

3.08 

1.55 

0.029 

0.88 

3.15 

1.62 

0.030 

0.90 

3.21 

1.68 

0.031 

0.92 

3.28 

1.75 

0.031 

0.94 

3.35 

1.82 

0.031 

0.96 

3.42 

1.89 

0.032 

0.98 

3.48 

1.95 

0.033 

1.00 

3.55 

2.02 

0.033 
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Since  one  should  be  integrating  for  the  optimal  time,  the  one  should  choose  should 
correspond  to  an  optimal  ,•  for  the  particular  i  h'  used.  Referring  to  Table  15,  such  a 
matching  pair  is  i  A'  =  0.031  and  ;•  =  0.92.  Consequently,  the  scanning  loss  for  61  pulses 
is  about  1.75  db.  This  procedure  was  repeated  for  various  v,  and  the  results  are  shown 
in  Table  18.  Hence,  both  the  optimal  integration  angle  and  scanning  loss  are  a  function 
of  the  number  of  pulses  integrated. 


Table  18 

Scanning  Loss  and  Optimal  Integration  .^gle  for  Various  \ 


Number 
of  Pulses 
Integrated 
.\ 

Scanning 

Loss 

(db) 

Optimal 

Integration 

Angle 

(radians) 

Number 
of  Fhjises 
Integrated 
.V 

Scanning 

Loss 

(db) 

Optimal 
Integration 
Angle  •, 
(radians) 

61 

1.75 

0.92 

201 

1.55 

0.91 

101 

1.63 

0.92 

401 

_ _ _ 

1.42 

0.88 

Optimal  Weighting  (U) 

In  the  previous  subsection,  the  optimal  integration  angle  for  a  uniform  weighting  was 
found  for  the  case  of  a  radar  scanning  past  a  point  target;  i.e.,  the  amplitudes  of  the  re¬ 
turning  pulse  train  are  proportional  to  the  square  of  the  one-way  antenna  voltage  gain 
pattern.  In  this  section,  the  optimal  weighting  function  is  numerically  calculated. 

To  find  the  optimal  weighting,  let  us  consider  the  following  binary  hypothesis: 

11^:  the  z  consist  of  noise  along,  :  o, 

//,:  the  .  consist  of  signal  plus  noise,  .s  -  ■> , 

where  {2,  >  are  samples  from  the  output  of  an  envelope  detector  and  v  is  the  amplitude  of 
a  sinusoidal  signal.  Rice  (12)  calculated  the  probability  densities  to  be 

/■ :  z  ^  01 

and 

/-■'Z, !  >■  /■  o  t 

The  Neyman -Pearson  Lemma  states  that  the  optimal  decision  rule  should  be  based 
on  the  likelihood  ratio  /,,/Z,  i  given  by 


z 


Z.  2. 


-■Zr-^ 


/.,/Z, 


77  r 


Z  (- 


1 

i 

j 


which  reduces  to 
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^  2  , 


Since  the  exponential  term  is  a  constant,  the  likelihood  ratio  is  a  monotonlc  function 
of  the  product  of  the  Bessel  functions: 


77 


(45) 


and  this  product  is  the  optimal  detector.  However,  this  optimal  detector  is  never  imple¬ 
mented  because  of  its  cumbersome  form.  Rather,  the  usual  approach  is  to  expand  the 
Bessel  function  in  an  infinite  series; 


/  'J  '  ::  1  T  -  •  - 

16  1 


Substituting  Eq.  (46)  into  Eq.  (45)  yields 


(46) 


•  77i'  ’7-^ 


which  can  be  rewritten  as 


I  -  V 


w  V 

1  V’  V 


Uftv! 


(47) 


The  usual  small-signal  argument  is  now  employed;  that  is,  if  the  signals  are  small, 
i.e,,  if  ’  is  of  the  order  of  ,  the  higher  order  terms  of  Eq.  (47)  can  be  ne¬ 

glected  and  only  the  single  summation  need  be  considered; 


L, V  /rv;  . 


(48) 


Equation  (48)  represents  a  fourth-power  voltage -antenna-pattern  weighting,  since  s, 
is  the  square  of  voltage  antenna  pattern.  This  is  in  contrast  to  the  usual  weighting  that 
is  used  — the  squared  weighting.  Fortunately,  much  of  the  difficulty  can  be  removed  if 
one  realizes  that  Eq.  (48)  is  not  a  very  good  approximation  to  Eq.  (47)  unless  the  signals 
are  extremely  small.  This  is  because  even  though  the  fourth-power  terms  are  of  the 
order  there  are  i/-  of  these  terms  in  comparison  to  the  '/  terms  of  order  .  Thus, 
the  optimal  weighting  ,u  for  a  square  law. 


(49) 


is  not  ,  since  the  fourth  power  terms  which  cannot  be  neglected  depend  on  the  squared 
terms.  Consequently,  it  will  be  necessary  to  calculate  . 


SECRET 


38 


G.  V.  TRUNK 


SECRET 


The  method  of  calculation  is  as  follows;  First,  if  expression  (49)  is  used  as  our  test 
statistic  and  if  the  calculations  in  the  previous  section  are  repeated,  the  probability  of 
detection  can  be  shown  to  be  a  monotonic  increasing  function  of  the  quantity 


/ 


:  =  - 


2  ~  -  <]>■  ‘  (  1  -  P,  , ) 


2 


=  -.\ 


(50) 


The  optimal  weighting  pattern  consists  of  the  set  of  ■  '  which  maximizes  Eq.  (50).  The 

optimal  weighting  was  found  by  numerical  methods,  and  the  results  are  shown  in  Fig.  13. 
As  one  can  readily  see,  the  optimal  weighting  is  neither  the  square-law  or  fourth-law  but 
rather  something  in  between.  Also,  the  weighting  is  a  function  of  the  n  v  ratio  per  pulse. 

The  question  naturally  arises  as  to  how  much  improvement  one  obtains  by  using  the 
optimal  weighting  instead  of  a  sin  weighting  or  a  uniform  weighting.  As  summa¬ 
rized  in  Table  19  the  results  indicate  that  if  one  integrates  over  an  interval  ,  the  opti¬ 
mal  weighting  is  only  0.1  db  better  than  the  uniform  weighting.  Also,  if  one  doubles  the 


Fig.  1  J 


\\'i  ighiii^g 


C'  I  lO  1.  s 
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integration  interval  (doubles  the  storage  re¬ 
quirements),  the  optimal  weighting  is  only  0.3 
db  better  than  the  uniform  weighting.  Thus, 
the  uniform  weighting  over  the  angle  ,  is 
very  close  to  the  optimal  condition. 


Azimuthal  Position  Estimates  (S) 

The  question  considered  in  this  subsection 
is  the  following:  After  a  target  has  been  de¬ 
tected,  how  close  can  its  position  be  estimated 
in  azimuth?  The  two  estimators  of  position 
Which  are  considered  are  shown  in  Fig.  14. 
The  first  estimator  is 

X  ^  .|\1A\, 

where  J\l\x  is  that  j  such  that 


Table  19 

Output  .«  v  Ratio  for  Different 
Weighting  Patterns 


Weighting 

Pattern 

■A  N  (db)  for 

Two  Lengths  of 
Integration 
Interval" 

Uniform 

15.7 

14.6 

<  s  i  n  J"  J- 

15.77 

15.9 

Optimal 

15.8 

16.0 

The  ii'.p-.it  .S  \  ratio  was  1.4  db  per 
pidse,  and  i;  tlie  opt  ima  1  inte  pr  a - 
tion  angle  for  a  iiniforrri  v.ei  ght  mg  . 


and  the  second  estimator  is 


-  1  .11  \si  r  JMItSil  J, 

where  Ji  ins  i  is  the  smallest  \  such  that  >  r  and  ji  \s  i  is  the  largest  j  such  that 

>’,./)  11. 


value  of  ) 

Fie.  14  _  A/.inmtha!  po^ilior.  .'stimates: 
)  \1  A  \  and  .  1 1  \>  I  •  .1 1  I  ll>ri  C 


While  it  is  quite  difficult  to  calculate  the  distributions  of  .■  and  '  ,  a  Monte  Carlo 
study  consisting  of  1000  cases  shows  that  '  and  are  both  gaussian  unbiased  estimators. 
This  inference  can  be  arrived  at  from  Figs.  15  and  16,  in  which  a  straight  line  indicates 
a  gaussian  distribution.  Hence  the  variances  oi  the  estimators  indicate  the  relative 
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Fi.’,  15  -  Dislribvuion  of  /  I  M  \\ 


goodness  of  the  estimators.  The  standard  deviations  tor  a  1000-case  Monte  Carlo  study 
are  shown  in  Table  20.  All  the  estimators  are  quite  compatible  in  pertormance  e.xcept 
for  ^  with  a  uniform  weighting.  One  can  make  the  following  statement  about  azimuthal 
position;  The  interval 

j-  -  J  (  \  .1  r  ■  ;  ‘  .  ■  !  t>  J  *  J  V  a  I  7  ‘ 

or  the  interval 

’■  -  \  ,1  r*  ;/»  I  '  -  II.  V  •  -  '  \  .1  r  ■  *;  '  I  '  “ 


Table  20 

Standard  Deviations  of  the  Estimators  and  v 
for  a  1000-Case  Monte  Carlo  Study 


Estimator 

Standard  Deviation  for  Different 
Weighting  Patterns 

_ 

Uniform 

,  1 

_ _ 

Optimal 

' 

0.20 

0.15  , 

0.15  , 

0.15  , 

0.13  , 

0.13  , 

■  ,  iM  iIk-  .tnulc 

to  tliC  :i.i If-powt  r  poir.l. 

SECRET 


SECRET 


NRL  REPORT  6804 


41 


Fig.  1 6  -  Distribution  of  J  -  t  JIIIISP  J 


covers  the  true  target  azimuthal  position  95  percent  of  the  time.  This  means  that  witli  a 
beamwidth  of  0.4  degree  ,  -  0.4  degree)  and  at  a  range  of  2000  naut  mi,  95  percent  of 
the  time  the  target  can  be  located  within  the  accuracies  given  in  Table  21.  The  general 
formula  for  accuracy  is 


2ii-:  >. 

where  R  is  the  range  in  nautical  miles  and  ■  ‘  "  ,1  is  the  standard  deviation  of  the  estima¬ 
tor  (given  in  Table  20)  in  terms  of  ,  in  radians,  with  ,  being  the  beamwidth. 


Table  21 

Azimuthal  .Position  Accuracy  With  Which  a  Target 
Can  be  Located  95  Percent  of  the  Time 


Azimuthal  Position  Accuracy  (naut  mi) 

Estimators 

for  Different  Weighting  Patterns 

Uniform 

1 s  1  n  /•  ’■  ./  J  ■ 

Optimal 

X 

2.79 

2.09 

2.09 

y 

2.09 

1.81 

1.81 
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All  of  the  previous  statements  on  azimuthal  accuracy  have  been  based  on  the  as¬ 
sumption  of  an  integrated  .'  \  latio  of  about  16  db.  Obviously,  larger  targets  will  have  a 
much  larger  integrated  A  ratio,  and  then  better  azimuthal  accuracies  result.  For  in¬ 
stance,  for  an  output  s  \  ratio  of  26  db  with  a  uniform  weighting  the  standard  deviations 
of  7  and  ,y  are 

V,,'.;'  -  -  0  oTi  .  (51a) 

and 

.\  ,1  r-'  .v  )  ;  ■  -  0  tU'-,  .  (51b) 

Using  Eqs.  (51),  the  resulting  azimuthal  accuracies  are  1.02  naut  mi  for  ?  and  0.C8  naut 
mi  for  v  for  an  0.4-degree  beamwidth  and  a  2000-naut  mi  range. 

The  question  arises:  What  is  the  maximum  accuracy  obtainable  for  the  previous  set 
of  conditions?  Fortunately,  tliis  question  can  be  answered  by  making  use  of  some  previ¬ 
ous  work  of  Swerling  (13).  He  showed  (by  using  the  Cramer  Rea  inequality)  that  the 
standard  deviations  of  the  optimal  estimator  were  ,  =  0.11  for  the  16-db  case  and 
.  ,  =0.031  for  the  26-db  case.  These  accuracies  correspond  to  position  errors  of  1.53 

naut  mi  for  16  db  and  0.43  naut  mi  for  26  db.  Thus,  the  two  simple  estimators  are  not 
too  tar  off  from  the  optimal  estimator. 


METHOD  OF  MAINTAINING  A  CONSTANT  FALSE  ALARM  RATE  (U) 

Consideration  of  Different  Methods  (U) 

It  is  usually  assumed  that  the  noise  has  a  Rayleigh  distribution.  Unfortunately,  in 
the  surveillance  problem,  the  radar  will  be  limited  by  sea  clutter  at  the  large  grazing 
angles  as  opposed  to  receiver  noise  (Rayleigh)  limited  at  the  shallow  grazing  angles. 

While  one  usually  assumes  a  Rayleigh  distribution  for  the  sea  clutter,  it  appears  very 
unlikely  that  tlie  noise  has  a  Rayleigh  distribution  out  to  6  .  (This  corresponds  roughly 
to  a  threshold  for  a  =  10" '  ' .)  Consequently,  to  maintain  a  /'  =  10' '  '  when  one  is 

clutter  limited,  one  usually  either  uses  aiaplive  schemes  or  else  hard  limits  the  data. 
Usually  adaptive  methods  are  used  w'hen  the  form  of  the  noise  distribution  is  known  but 
certain  of  its  parameters  are  unknowm.  For  instance,  if  the  noise  is  gaussian  with  a 
nOnstationary  unknown  variance  (noise  power),  the  adaptive  procedure  is  nothing  more 
than  a  variable  threshold  which  is  proportional  to  the  estimated  noise  power.  Since  the 
distribution  ol  sea  clutter  is  not  only  non-Rayleigh  but  also  unknown,  before  conventional 
adaptive  methods  could  be  used,  one  would  have  to  estimate  the  functional  form  of  the 
density  of  sea  clutter.  Estimating  this  density  out  to  6  or  <  seems  impossible.  Hence, 
adaptive  methods  of  this  nature  do  not  seem  very  promising. 

On  the  otlier  liand,  liard  limiting  seems  quite  practical  il  one  is  integrating  for  .a 
longtime,  i.e,,  accumulating  25  or  more  pulses.  While  liard  limiting  yields  a  constant 
false  alarm  rate,  it  results  in  a  loss  ol  about  2  db.  (Further  discussion  ot  the  hard  limit¬ 
ing  loss  is  given  in  Ret.  14.)  Untortunately,  for  certain  types  of  processing,  hard  limiting 
also  causes  strong  targets  to  interfere  with  the  detection  of  weaker  targets  (15). 

Thus,  both  adaptive  and  liard  limiting  methods  have  undesirable  features.  However, 
a  method  comihning  soft  limiting  and  an  adaptive  threshold  seems  quite  promising.  Thi.s 
method  is  based  on  the  fact  tliat  even  though  the  density  of  sea  clutter  is  nongaussian  in 
the  6  or  7  region,  it  is  almost  gaussian  below  2.5  or  3  .  Hence,  if  one  limits  the  data 
at  3  (for  instance),  one  now  knows  the  density  of  sea  clutter  except  for  the  noise  power 
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-  .  Thus,  one  can  estimate  and  use  this  estimated  \alue  to  control  the  value  of  the 
soft  limiter  and  the  decision  threshold.  An  alternate  way  of  accomplishing  the  same 
thing  IS  to  have  a  fixed  value  for  the  soft  limiter  and  to  have  a  fixed  threshold.  The  es¬ 
timated  is  used  to  control  an  age  so  that  the  noise  is  always  limited  at  its  3  point.  An 
added  advantage  of  a  soft  limiter  is  that  it  reduces  the  dynamic  range  required  by  the 
data  processor. 


Soft  Limiting  Loss  (U) 

While  the  hard  limiting  loss  is  fairly  easy  to  calculate,  tlie  soft  limiting  loss  is 
much  more  difficult  to  calculate.  Marcum  (16)  defines  the  soft  limiting  ratio  as  the 
ratio  of  the  limit  level  to  the  rnis  noise  level.  He  concludes  that  if  the  number  of  pulses 
integrated  is  large,  the  limiting  loss  is  only  a  fraction  of  a  db  if  the  limiting  ratio  is  as 
large  as  2  or  3,  but  that  if  only  one  or  two  pulses  are  integrated,  the  limiting  ratio  must 
be  in  tlie  neighborhood  of  10  to  prevent  a  serious  loss.  Thus,  while  Marcum  makes  the 
qualitative  statement  that  in  certain  cases  the  loss  is  small,  he  makes  no  quantitative 
statement  about  the  loss,  the  limiting  value,  or  the  number  of  pulses  integrated.  Conse¬ 
quently.  it  was  decided  to  investigate  tliis  problem  further. 

In  attacking  tliis  problem,  one  could  use  the  standard  approach  that  employs  the 
Gram-Charlier  series;  however,  Marcum  said  the  calculations  were  quite  tedious.  Al¬ 
ternately,  one  could  calculate  the  probability  distribution  of  the  limited  noise  by  the 
characteristic  function  method.  Unforlunaiely.  this  results  in  a  triple  integral  for  which 
there  c-.xists  no  known  closed  lorni  solution.  A  numerical  integration  of  this  triple  inte¬ 
gral  to  the  desired  accuracy  would  require  an  cxoi'bitaiit  amount  of  lime  on  a  digital 
computer;  consc  'iently,  this  method  seems  impractical  also.  Therefore,  it  was  decided 
to  use  a  Monte  .  lo  method  to  calculate  the  loss  rather  tlian  use  cither  of  the  two  pre¬ 
vious  mctl'.ods. 

The  system  to  be  investigated  by  the  Monte  Carlo  (iroccdure  is  shown  in  Fig.  17.  It 
consists  oi  a  .soft  limiter,  followed  by  a  dechirping  device,  followed  liy  an  envelope  de¬ 
tector.  lollowed  Ijv  a  b-pulse  accumulator  wliose  output  is  cumpaixd  to  a  threshold  ’  . 
The  Monte  Carlo  i/i  ocedurc  just  consists  in  iierlorming  the  operations  shown  in  Fig.  18 
wliich  IS  tlie  mailiemalical  model  of  tlie  system  shown  in  Fig.  17.  The  main  part  of  the 
simulation  is  (lie  gcnei'iition  of  the  signal  plus  noise  components  ,  and  .  '  .  The 
gaus.«i.Ui  noise  .  is  generated  by  the  equation 


where  .  is  Rayleigii  distribuiod  and  is  uniformly  distribuicd  on  (0.  2  ).  Since  a  uniform 
random  number  generator  is  available  on  the  computer  wiiu  li  geiicraie.s  a  random  num¬ 
ber  uniformly  distributed  on  (0,  1),  generation  of  simply  involves  the  multi  plication 
of  .  by  '  .  The  goneratioii  of  is  accomi  hslied  by  the  uniform  random  number  gener- 
<itoi'  in  liif  lollowiiig  way: 
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Fig.  18  -  MathcDKit  ica  1  niorli  l  of  lla- 


Integrating,  one  obtains 


or 


/'(  /•  ./) 


I  -  !•< 


\  ■>.- 


The  quantity  l  -  /’(  /■-  is  uniformly  distributed  between  :  and  1  and  ;  oiisociueni ty  can 
be  replaced  by  a  random  numlter  r .  Hence,  if  one  S()lve.‘^  llq.  (52)  for  .  one  obtains 


,1  (  -C  1m  I  )  '  ' 

which  shows  that  n  has  a  Rayleigh  distribution.  Thus,  one  can  gonerati  the  component.^ 
j  '.  -  and  y'.  by  generating  two  random  numbers  -  .  and  •  ;  i.e., 


and 


where 


and 


J-  -  .  .S  I 


V  •  •  f/  ■  «  ' 


i".  .  2 

'  } 


"i  j  '  -  I',  r,  , 


To  find  the  limited  components  x- j  and  one  first  calcuhdcs  the  envelope  >  : 
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If  c ,  ,  1  SI.  (soft  limiting  value),  then  j  .  j,'  .  and  v; ,  =  y,'. .  On  the  other  hand,  if 

f,.  >  SI  ,  then  X.  .  =  SI  (2-',)  c,.  and  y, .  -  si. . . 

After  the  soft  limiter,  the  simulation  process  is  straightforward.  One  generates  \ 
pairs  of  pulses  in  a  similar  manner,  sums  up  these  v  pulses,  squares  the  two  sums,  adds 
the  two  squares,  and  takes  the  square  root  —  this  represents  the  output  of  an  envelope 
detector  of  a  dechirped  pulse.  Finally,  one  sums  i/  of  the  dechirped  pulses  and  compares 
the  sum  to  a  threshold  l . 

Finally,  the  system  was  simulated  on  the  CDC  3800  with  the  following  input  parame¬ 
ters:  \  =  63.  .'/  =  128,  =  0.13515,  =  1,  and  T  =  1620.  The  threshold  was  set  to  yield  a 

probability  of  false  alarm  of  lO'"  and  a  probability  of  detection  of  0.5,  The  results  of 
the  Monte  Carlo  procedure  are  shown  in  Table  22. 


Table  22 

Results  of  the  Monte  Carlo  Procedure  Shown  in  Fig.  18 


Limit 

Value 

No.  of 
Cases 

No  Of 

Cor  .-ect 
Decisions 

Probability 
of  Detection 

Estimated 
t  Ratio 

1 

1 

Loss 

(db) 

107 

57 

.  '  ■  ' 

0.532 

-2.3 

3.0 

107 

49 

0.458 

-2.5 

0.2 

2.5 

107 

27 

0.252 

-3.0 

0.7 

2.0 

59 

0.018  1  -4.0 

1.7 

1 

The  100  cases  or  so  that  were  run  provide  a  good  estimate  of  the  probability  of  cor¬ 
rect  decision.  However,  one  billion  cases  would  have  to  be  run  if  one  wanted  to  calculate 
the  probability  of  false  alarm.  But,  since  the  inequality  T.  ■  10"'  is  true,  one  can 
bound  the  soft  limiter  loss  by  setting  ,  =  10  "  for  all  limit  values.  The  estimated  \ 
ratios  are  found  by  finding  the  s  \  ratios  tlial  yield  the  i\  values  in  Table  22  when  one 
integrates  128  pulses.  The  v  \  ratios  were  found  using  Robertson's  curves  (17).  Con¬ 
sequently,  the  losses  for  limiting  values  of  3  ,  2.5  ,  and  2  are  bounded  by  0.2,  0.7,  and 
1.7  db  respectively.  I'hus,  as  long  as  one  is  limiting  above  2.5  ,  the  loss  is  lower  than 
0.7  db;  and  it  is  this  author's  opinion  that  the  actual  loss  is  abo".l  0.2  db. 
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Appendix  A  (U) 

DERIvyiTTON  OF  EQ.  (42) 

Given  two  gaussian  distributions,  (i~  >  and  •.■i  ,  representing  the  null 

hypothesis  and  the  alternative,  the  probability  of  false  alarm  is 


and  the  probability  of  detection  is 


(Al) 


-  j  '.'  r  ,  (A2) 

where  'I  is  the  threshold.  Equations  (Al)  and  (A2)  can  be  rewritten  as 

o'to,  i  'hi'  7  -  •-■I'i  (A3) 

O'-.  o,n.;v.  (A4) 

.  t  7 

One  may  imply  from  Eq.  (A4)  that  is  a  monotonic  increasing  function  of  the  quantity 

7';  =  7'  7't.  (A5) 

Solving  Eq.  (A3)  for  / ,  one  obtains 

I'  -  •  |'1‘  -  I  -  i\  (A6) 

Substituting  Eq.  (A6)  into  Eq.  (A5)  yields 

I',.-  --  7  I  i...  -  -  ,  'I'".  I  -  7’,^.)  I  ,  (A7) 

where  7'  is  a  monotonic  increasing  function. 


and 


I  -  , 


7’; 


'f 
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Several  statistical  topics  were  studied  concerning  detection  of  targets  on  the 
ocean's  surface  from  a  satellite-based  radar.  First,  a  desired  signal-to-noise 
ratio  was  obtained  by  considering  the  tradeoff  between  false  alarms  and  detected 
targets.  For  a  0.99  detection  probability  and  a  10"  false  alarm  probability  the 
S  N  ratio  would  be  16  db  for  a  nonfluctuating  target.  Second,  the  decorrelation 
times  of  sea  clutter  were  calculated  and  found  to  be  so  small  that  for  side-looking 
radars  sea  clutter  samples  are  independent  from  pulse  to  pulse  but  found  to  be  so 
large  for  forward-scanning  radars  that  one  finally  obtains  decorrelation  from  sea 
motion  instead  of  from  platform  motion.  Third,  a  fluctuating  target  was  considered, 
and  the  losses  were  calculated  and  compared  with  those  of  a  nonfluctuating  target. 
Roughly,  one  needs  1  to  4  db  more  power  to  detect  a  fluctuating  target  as  opposed  to 
a  nonfluctuating  target.  Then,  the  ionosphere  was  considered,  and  the  effects  of 
Faraday  rotation  and  the  random  phase  shift  w-ere  calculated.  Degradation  varies 
from  intolerable  below  900  MHz  to  negligible  at  3  GHz.  Next,  some  of  the  aspects 
of  data  processing  were  studied;  the  optimal  integration  angle  was  found,  the  opti¬ 
mal  weighting  was  calculated,  and  azimuthal  position  estimators  were  considered. 
Surprisingly,  optimal  weighting  gives  a  pattern  only  0.3  db  better  than  uniform 
weighting.  Finally,  a  way  was  determined  of  using  an  adaptive  threshold  and  soft 
limiting  to  maintain  a  constant  false  alarm  rate. 
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